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ABSTRACT 

Before reionization, the intergalactic medium (IGM) may have been sufficiently cold for low-mass 
"minihalos" to condense out of the gas and subsequently affect reionization. Previous work has shown 
that minihalos generate reasonably large 21 cm fluctuations. Here we consider this signal in its 
proper cosmological context and show that isolating minihalos from the rest of the IGM is extremely 
difficult. Using the well-known halo model, we compute the power spectrum of 21 cm fluctuations 
from minihalos and show that the signal decreases rapidly as feedback increases the Jeans mass. We 
then show that even a small Lya background increases the 21 cm fluctuations of the diffuse IGM well 
beyond those of the minihalos; because the mass fraction in the IGM is much larger, minihalos will 
lie buried within the IGM signal. The distinctive signatures of non-linear bias and minihalo structure 
emerge only at much smaller scales, well beyond the resolution of any upcoming instruments. Using 
simple, but representative, reionization histories, we then show that the required Lya background 
level is most likely achieved at z > 15, while minihalos are still rare, so that they are almost always 
degenerate with the diffuse IGM. 

Subject headings: cosmology: theory - intergalactic medium - diffuse radiation 



1. INTRODUCTION 

Cosmologists have long recognized that the ubiq- 
uity of hydrogen in the intergalactic medium (IGM) 
makes its line transitions particularly effective diag- 
nostics of the diffuse material that separates galaxies. 
One especially useful prob e is the so-called Lya for- 
est IjGunn fc Petersonl H965'). which has proved exceed- 
ingly useful for studying the z < 6 Universe. Un- 
fortunately, this transition is so strong that the for- 
est saturates at relatively small neutral fractions (e.g., 
iFan et all 12005(1 . Thus the much weaker 21 cm hy- 
perfinc transition is the probe of choice for study- 
ing the predominantly neutral gas in the high-redshift 
IGM, especially during and before reionization (Field 
1959 al ISunvaev fc Zcldovic hl Il975t iHoean fc Reeslll979t 
IScott fc Tta*«199ttlM«faii et al.Hl 997ft. 

Because it is a spectral line, 21 cm observations offer 
the possibility of mapping the three dimensional distribu- 
tion of gas from z ~ 200, when the gas temperature first 
decoupled from the CMB, to z ~ 6, a period spanning the 
"Dark Ages," the onset of structure formation, the birth 
of the first galaxies, and reionization itself. One particu- 
larly interesting application is to study the first collapsed 
objects. The IGM is visible only so long as the 21 cm 
spin temperature T$ differs from the CMB temperature 
T 7 (because the latter is used as a backlight). At z < 50 
(and before luminous sources appear), the vast majority 
of gas is invisible because collisions are not rapid enough 
to drive T$ — * Tr-, the kinetic temperature of the gas 
(which is colder than T 7 because of adi abatic expansion) . 
Only overdense regions remain visible. Illiev et aLl (2002, 
2003) pointed out that "minihalos," which are dark mat- 
ter halos larger than the Jeans mass but too small to 
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cool and form stars, constitute a particularly interesting 
class of 21 cm emitters. These objects are a qualita- 
tively new feature of the pre-reionization IGM that can 
slow r eionization because of their rapid recombination 
rates llHaiman et alJl200H lihev et alJl2003 iCiardi et all 



Although they are much too small to be imaged 
individually in realistic experiments, they imprint large- 
scale fluctuations on the 21 cm background that could, in 
principle, be used to c onstrain cosmological parameters 
(Illiev et alJl2lMl200l . 

The purpose of this paper is to consider the mini- 
halo signal in more detail and to embed it in its proper 
cosmological context. Because minihalos cannot be re- 
solved, they are obvio usly subject to confusion w ith other 
sources in the beam. IFurlanetto fc Loea (|2004f) pointed 
out one possibility: collisions can be sufficiently rapid in 
moderately overdense IGM gas to render visible sheets 
and filaments in the cosmic web. This has since been 
confirmed in numer ical simulations l|Kuhlen et al.ll200a: 
iShapiro et alJ l2005|) , where sheets and filaments con- 
tribute a significant fraction of the signal at z > 15 
and confuse attempts to identify the minihalos. But 
here we will focus on the low-density, diffuse IGM. 
The Wouthuysen-Field mechanism, in which absorption 
and re-emission of Lya photons mixes the spin states 
l|Wouthuvsenl 119521: iFieldl 11953 Iffiratal 120051 can drive 
Ts —> Tk everywhere and hence render the entire IGM 
visible. IQh fc Mackl l|2003j) argued that, if Wouthuysen- 
Field coupling is strong, the minihalos (which contain 
only a small fraction of the mass) will be buried in the 
rest of the signal. 

A related question is how efficiently minihalos can 
form once feedback from the first luminous sources is 
taken into account. In particular, X-ray heating in- 
creases the cosmological Jeans mass. The resulting 
"entropy floor" prev ents gas accretion onto minihalos 
(|Oh fc Haimanl 120031 hereafter OH03), effectively sup- 
pressing their abundance. But the importance of this 
mechanism is unclear, because X-rays can also stimu- 



late cooling, and entropy injection will be less efficient 
in minihalos that have already begun to collapse. One 
of the key questions of reionization is the efficiency of 
feedback, and 21 cm observations of minihalos could po- 
tentially shed light on it. Here we will develop a new 
method to calculate the 21 cm fluctuations of miniha- 
los (based on the halo model) and use it to predict the 
observable implications of feedback. 

The rest of this paper is organized as follows. We be- 
gin in i]21 and JU with brief discussions of minihalos and 
21 cm physics. In 21 we compute the fluctuation ampli- 
tude expected from minihalos as a function of the IGM 
entropy. We then relate that signal to fluctuations from 
the diffuse IGM in ^SJ including Lya coupling. Next we 
consider minihalos in the context of some representative 
reionization histories in Sj^l Finally, we conclude in jJJ| 

In our numerical calculations, we assume a cosmol- 
ogy with n m = o.26, n A = o.74, n b = 0.044, h = 

lOO/i km s" 1 Mpc" 1 (with h = 0.74), n = 0.95, and 
as = 0.8, c onsistent with the most recent measurements 
IjSpergel et al.ll2006j) . 3 Unless otherwise specified, we use 
comoving units for all distances. 

2. MINIHALOS AND FEEDBACK 

While collisionless cold dark matter can collapse into 
gravitationally bound objects of arbitrary masses, ther- 
mal pressure prevents baryons from accreting onto a halo 
unless its mass exc eeds the time-average d Jeans mass, or 
the "filter" mass l|Gnedin fc Huilll998D . We take this, 
evaluated in an adiabatically cooling universe, to be the 
minimum allowable mass m m i n for a minihalo. We set the 
maximum mass m max to corresp ond to a virial tempera - 
ture T vir = 10 4 K (defined as in lBarkana fc Loebll200ij) . 
the point at which atomic cooling becomes efficient. 

Of course, if the IGM temperature exceeds its nominal 
value from adiabatic cooling, the Jeans mass increases 
and minihalo formation is suppressed. OH03 showed that 
the IGM "entropy" ^igm = T K /n 2 / 3 quantifies the sup- 
pression; here Tr- is the kinetic temperature of the IGM. 
Clearly K is conserved so long as the gas evolves adia- 
batically which occurs both during the overall expansion 
of the universe and during the initial stages of halo col- 
lapse. Subsequent shock heating can further raise the 
gas entropy, but so long as radiative cooling is unimpor- 
tant, K\qm provides a well-defined entropy floor. We 
can similarly define K m h(r) = T vir /n(r) 2 / 3 to be the en- 
tropy profile acquired by a minihalo during gravitational 
collapse. We assume that the gas densi ty profile n(r) 
traces the universal dark matter profi le of iNavarro et al.l 
(J1997D . Numerical simulations fe.g.. JFrenk et al.lll999^ 
show this to be a good assumption except near the halo 
center, where in any case the effects of a finite entropy 
floor predominate. 

With this entropy profile, it is possible to calculate 
detailed equilibrium density and temperature profiles. 
Indeed, by constructing these profiles as a function of 
Kigm, OH03 showed that the structure and total ac- 
creted gas fraction f g (relative to the cosmic mean 

Qb/^m) are simply functions of K = Kigm / K m h(r v i T ) 
(see their Fig. 6). If K <C 1, the entropy generated by 
the gravitational accretion shock vastly exceeds the IGM 

3 Note that we have increased as above the best fit WMAP value 
in order to improve agreement with weak lensing measurements. 



entropy and f g — * 1. If, on the other hand, K > 1, the 
thermal pressure of the IGM exceeds the minihalo's grav- 
itational attraction, so f g — > 0. We compute f g for each 
minihalo by interpolating the results of OH03. N ote that 
these gas fractions match the fitting formula of iGnedml 
( 2000) reasonably well; the widely-used filter mass is sim- 
ply a crude estimate of the point at which K = 1 in an 
adiabatically cooling universe. 

This prescription for calculating density profiles (and 
hence f g ) assumes that K\qm is constant throughout the 
accretion history of the halo; in reality, i^iGM depends 
on the thermal history of the minihalo's progenitors, but 
we shall ignore such complications here. Furthermore (as 
in this paper), radiative cooling by H ? mole cules was as- 
sumed to be negligible. iRozas et alJ l)2005|) obtain com- 
parable results with a similar model. Note that if the 
soft-UV background flux is small, the effects of radiative 
cooling, which erases the entropy floor established by an 
X-ray background, cannot be ignored ( Kuhlcn & Madau 
2005). Generally, such considerations only apply to ex- 
tremely hard spectra. 

3. THE 21 CM TRANSITION 

We review the relevant characteristics of the 21 
cm transition here; w e refer the interested reader to 
iFurlanetto et al.l IJ2006J) for a more comprehensive dis- 
cussion. The 21 cm brightness temperature (relative to 
the CMB) of a patch of the IGM is 
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where xhi is the neutral fraction, S is the fractional over- 
density, T 1 is the CMB temperature at redshift z, and 
Ts is the spin temperature. Note that the patch will ap- 
pear in absorption if Ts < T^ and emission otherwise. 
Here we have assumed that the patch expands uniformly 
with the Hubble flow; radial peculiar velocities also affect 
STb by changing the mapping from distance to frequency 
l|Bnaradwai fc Alil f2004: B arkana fc Loebll2005"al) . How- 
ever, these anisotropic fluctuations do not affect our ar- 
gument, so we will ignore them for simplicity. 4 

The spin temperature Ts is determined by the com- 
petition between three processes: scattering of CMB 
ffiotons, col lisions, and scat tering of Lya photons 
1 IT952llFiddi ri958V In equilibrium, 
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Here x a oc n is the collisional coupling coefficient 
IjZvgelmanl I2005J) . T K is the gas kinetic temperature, 
and T c is the effective color temperature of the radi- 
ation field. Normally f c w T K (|Fiek1 Il959bt IHiratal 
I2005|) . The middle term describes the Wouthuysen-Field 
effect, in which absorption and re-emission of Lya pho- 
tons mixes the hyperfine states. The coupling coefficient 

4 Note that the velocity field around an individual minihalo is 
far from the Hubble flow, so this is not a good assumption if the in- 
strument is able to resolve the halos. However, such high-resolution 
experiments will be impractical for the foreseeable future. 
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where S a is a factor of order unity describing the de- 
tailed atomic physics of the scattering process and J a 
is the background flux at the Lya frequency in units of 
cm~ 2 s _1 Hz -1 sr _1 ; the Wouthuysen-Field effect be- 
comes efficient when there is about one L ya photon per 
ten baryons. We use the numerical fits of Hirata (2005) 
for S a and T c . 

4. FLUCTUATIONS FROM MINIHALOS 

We will begin by calculating the power spectrum of 
the 21 cm brightness temperature assuming that only 
minihalos are allowed to emit. As we will see below, 
this requires an unlikely set of circumstances, but it pro- 
vides some intuition about the minihalo signal and the 
effect s of feedback upon it. We will use the halo model 
fe.g.. lCoorav fc Shcth 2002), which constructs the power 
spectrum of density fluctuations by dividing the universe 
into dark matter halos of all sizes. The power spec- 
trum has two components: a "one-halo" term that de- 
scribes correlations between particles within the same 
virialized halo and a "two-halo" term that describes cor- 
relations between particles in separate halos. For our 
purposes, the advantage of this approach is that the 
contribution from minihalos is obvious: we simply in- 
clude only those particles that lie within minihalos. Thus 
Pmh(fc) = PH(k) + P^(fc), where 
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Here 6Tb is evaluated via equation Q with Tg 3> T 7 and 
S = 0, p is the mean d ensity of the universe, rih (m) is the 
halo number density IjPress fc Schechter 1974), 5 u(k\m) 
is the Fourier transform of the normalized halo d ensity 
profile, b(m) is the linear bias l)Mo fc Whitdll996|) . and 
T is a temperature factor. Although the baryonic distri- 
bution differs slightly from the dark matter distribution 
because of its finite pressure, it suffices to let u(k\m) 
equal the lNavarro et al.1 l)1997l) profile for our purposes. 
As we shall see, this only affects the power spectrum on 
extremely small scales. 

The remaining factor is 7^, which is the temper- 
ature factor (Tc; — T^)jTs av eraged over each mini- 
halo. Illiev et all l|2002L l2003[) computed the detailed 
brightness temperature profiles of minihalos (see also 
iFurlanetto fc Loebl 120021) . They showed that, although 
Ts does vary across the minihalo, most gas still has 
T$ 3> T!y because it is so dense everywhere. In this case, 
Tmh ~ 1 • For simplicity, we will therefore evaluate "T^h at 
the mean overdensity of the minihalo (/ g A v ; r , where the 




5 Using the Shcth fc Tormeri 119991) mass function does not af- 
fect our conclusions. 



Fig. 1. — Expected rms brightness temperature fluctuations from 
minihalos. Thick solid and thin dashed curves are for z = 10 
and 20, respectively. From top to bottom within each set, the 
curves assume Tx = T a( j, 20, 100, and 1000 K (we set J a = 
throughout). The dot-dashed curves show the fluctuations from 
linear theory, assuming that the entire IGM has Tg ^> T 7 . 



latte r factor is the mean ov erdensity of a virialized ob- 
ject; |B^^^^^No™a^998). Note that this prescription 
includes (approximately) the effect of the entropy floor 
on the 21 cm signal; OH03 make more detailed estimates 
of the spin temperature profiles of individual minihalos 
by including the full equilibrium radial density profiles. 
This approximation does not affect our conclusions. 

Figure ^ shows the resulting rms temperature fluc- 
tuations as a function of scale, defined via A 2 (A:) = 
k 3 P n ih(k)/2ir 2 . The thick solid and thin dashed curves 
show results for z — 10 and z — 20, respectively. We as- 
sume J a = here, so that Lya coupling is insignificant. 
For reference, the upper thick and lower thin dot-dashed 
curves show the (linear) amplitude if the entire IGM had 
T s » T 7 . 

First consider the uppermost minihalo curve in each 
set. These assume that Tjc — T a( j, the temperature ap- 
propriate for adiabatic co oling (T„.^ = 2.58[( l + z)/ 10l 2 K 
according to RECFAST: (Seager et all Il999l 12000 K Our 
choice for m m ; n already includes the effects of the en- 
tropy floor for this temperature history, so f g ~ 1 in all 
minihalos. The signal has two components. At scales 
fc > 20 Mpc -1 , A increases rapidly. This comes from 
-fmh' an d ^ corresponds to scales on which the internal 
structure of minihalos is resolved. In reality, finite gas 
pressure will make this part shallower, but such scales 
are well beyond th e capabilities of any currently planned 
exper iments (e.g., IBowman et al.1 120051 iMcQuin n et alJ 
2005). Note tha t our version of th e halo model includes 
only linear bias. Illiev et al. ( 2003) have shown that non- 
linear bias modifies the power spectrum at fc > 6 Mpc -1 
- still somewhat beyond the capabilities of any of these 
observatories. On larger scales, the fluctuations trace the 
linear power spectrum, with the input parameters affect- 
ing only the amplitude. This is the regime where P^ 
dominates. On such scales, u(k\m) rs 1, so 7 m h — * 6/mh, 



where b is the average (mass- weighted) bias of the mini- 
halos and / m h is the fraction of the IGM mass incor- 
porated in minihalos. We find that (/ m h, b, I m h) ~ 
(0.067, 2.8, 0.19) and (6.7 x 10" 4 , 8.0, 0.0054) at z = 10 
and 20, respectively. Note that, because they are far out 
on the nonlinear mass tail, the minihalo bias evolves by a 
large factor over these redshifts. Nevertheless, minihalo 
fluctuations are considerably smaller than those from a 
uniformly hot IGM. We predict rms fluctuations < 2 mK 
(or 0.01 mK) at z = 10 (or 20) on scales accessible to ob- 
servations. The signal increases rapidly with cosmic time 
until z ~ 8, when minihalos become nearly linear fluctu- 
ations; at significantly higher redshifts they contain such 
a small fraction of the baryons that detecting them will 
be a true challenge 

These uppermost curves ignore the possibility of ther- 
mal feedback in the IGM. The other curves show the ef- 
fect of this "entropy floor": they take Tk = 20, 100, and 
1000 K, from top to bottom. At the minimum minihalo 
mass, these have K = (4.9, 24, 244) and (2.3, 11, 114) 
for z — 10 and z = 20. Obviously even a modest amount 
of heating dramatically reduces the minihalo signal, es- 
pecially at z = 10. This is because most of the mass 
is contained in halos with m « m. m in ; which are most 
sensitive to feedback. Thus the minihalo 21 cm signal 
is probably only detectable if the IGM is relatively cold 
at the time of their formation. Their sensitivity to Tk 
suggests that constraining the minihalo signal could be 
a powerful probe of feedback, but it also compromises 
attempts to use it to constrain cosmological parameters 
(jlliev et alJl2002j) . 

5. MINIHALOS AND THE IGM 

Of course, short of resolving minihalos, any observed 
patch of sky must contain a mix of low-density IGM gas 
and minihalos. To better mimic the actual observables, 
we now compute the fluctuations from minihalos together 
with those of the surrounding diffuse IGM. In the halo 
model, this IGM gas is made up of all halos with masses 
between zero and m m ; n . Because such halos are too small 
to accrete gas, the baryons remain at a constant density, 
so u(k\m) = 1. Additionally, the one-halo term vanishes. 
The temperature factor 7igm can be evaluated at the 
mean cosmic density, 6 given a temperature history and 
a Lya background. The total power spectrum of 21 cm 
fluctuations is then 

P 2 i(k) = P^(k) + 6T 2 b P lin {k)[I mh (k) + h GM (k)] 2 , (7) 

where 

IiGM(k) = Tigm 
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The second part is the usual two-halo term split into 
the IGM, minihalos, and the correlations between them. 
The temperature factors are crucial, because if Tk < 
T 7 the IGM absorbs CMB photons while the minihalos 
emit. In this case the two tend to cancel each other out. 
One consequence of equation J7|) is that, unless one can 
reach scales where P^ is large, the IGM and minihalo 

6 Actually, in the presence of an entropy floor, the IGM density 
is slightly larger than the cosmic mean because of gas that could 
not accrete onto minihalos; we assume that this gas is distributed 
uniformly throughout the universe. 




Fig. 2. — Expected rms brightness temperature fluctuations at 
z = 10 for J a = 0. The thin and thick curves have Tk = T a( j 
and Tk = 100 K. Within each set, the solid and dotted lines show 
the fluctuations from minihalos and the IGM, respectively, while 
dashed lines show the net result (for Tk = T a< j, this overlaps the 
solid curve). The dot-dashed curve shows the fluctuations accord- 
ing to linear theory, assuming that the entire IGM has Ts 3> T 7 . 



fluctuations are completely degenerate. Because both 
trace Pn n , any minihalo signal can be mimicked by a 
suitable choice for 7igm ■ Note that we have ignored halos 
with m > m max (which are able to cool rapidly and form 
stars). If they contain neutral gas, they would boost the 
signal slightly. 

Of course, our assumption of a constant density 
isothermal IGM is much too simplistic. Structure for- 
mation organizes the IGM into hot, overdense sheets 
and filaments that exhibit stronger coupling between the 
spin and kinetic t emperatures. Using a sim ple model of 
spherical collapse. iFurlanetto fc L ocb (2004) argued that 
when minihalos are still on the nonlinear tail of the mass 
function (at z > 15), these shocks will induce 21 cm 
fluctuations comparable to those of minihalos. Numer- 
ical simulations clearly show that these networks of fil- 
amentary emission produce interesting fluctuations that 
interfere with the minihalo signal, espe cially at z > 20 
(Kuhl en et al.l l2006: Sh apiro et alJl2005j) . and some show 
that Ts flu ctuations in the low-den sity gas can also be 
significant iJGnedin fc Shaver] I2004J) . Our purpose is to 
show that minihalos will be extremely difficult to sep- 
arate from fluctuations in the diffuse IGM; by neglect- 
ing its spin temperature and density fluctuations, we are 
therefore taking the most conservative possible stance. 

Figure compares the minihalo and IGM 21 cm fluc- 
tuations at z = 10, assuming J a — 0. First consider the 
thin curves, which assume Tk = T a< \. The dot-dashed 
curve shows the linear density fluctuations if Ts S> T r 
The solid curve shows the minihalo component, and 
the dotted curve shows the IGM component. Minihalos 
clearly dominate the signal: the IGM temperature and 
density are such that x c <C 1 at z — 10. However, if we 
increase the IGM temperature to Tk = 100 K, as shown 
in the thick curves (here the dashed line shows the net 
fluctuations), the situation changes. First, the minihalo 



fluctuations decrease by an order of magnitude because 
the entropy floor suppresses accretion onto them. At 
the same time, the IGM fluctuations strengthen because 
x c increases with temperature, even though the density 
remains (nearly) the same. Thus if Tk > 100 K, the 
observed fluctuations no longer directly trace the mini- 
halo population, even if J a = (provided feedback is as 
strong as our model predicts). At higher redshifts, the 
crossover occurs with even less heat input because the 
fraction of gas in minihalos is smaller. Clearly the possi- 
bility of IGM heating can significantly reduce our power 
to interpret weak 21 cm fluctuations as minihalos. 

Thus far we have assumed J a = and included only 
collisional coupling. Figure [3] shows how a Lya back- 
ground affects the IGM signal at z = 10 in a sce- 
nario with no X-ray heating. We emphasize that this 
maximizes the minihalo component by ignoring the en- 
tropy floor. The solid curve shows the minihalo con- 
tribution (which is essentially independent of J Q , be- 
cause collisional coupling is always efficient at virial den- 
sities and temperatures). The dashed curves show the 
net signal (including the IGM) assuming J a = 2 x 
(10~ 12 , 10~ h , 10- 10 ) cm" 2 s- 1 Hz" 1 sr- 1 (x a ~ 0.04- 
4), from bottom to top. For contex t, the ionized fraction 
at a given x a is l|Furlanettoll200 6b) 
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,l + n rec 0.1 N a S a J\l + z 

where / osc is the escape fraction of ionizing photons, 
ri rec is the mean number of recombinations per hydro- 
gen atom, and N{ on /N a is the mean number of ionizing 
photons per Lya coupling photon. This is rs 0.45 for 
Population II stars and rs 7 for very massive Population 
III stars. We thus expect Lya coupling to turn on well 
before reionization is complete. Also note that a soft-UV 
background of x a ~ 0.002-0.02 su ffices to suppress H 2 
formation and cooling in minihalos (Haiman et al. 2000), 
so in principle there can exist a period during which mini- 
halos cannot cool but dominate the 21 cm emission. 

Figure |21 shows that, once Lya coupling begins, the 
IGM rapidly overwhelms the minihalo signal. This is 
because the IGM contains so much more mass than the 
minihalos and because (unlike emission) absorption does 
not saturate l)Oh fcMackl 12003) . We find that, once 
x a > 0.02, the net fluctuations no longer closely trace 
those of minihalos for x a < 0.002, the minihalos domi- 
nate. As suggested by equation 10 and confirmed in the 
next section, this occurs well before reionization. Thus 
minihalos will only be visible at quite high redshifts. The 
Lya background is the most likely mechanism to render 
minihalos effectively invisible. Note as well that the IGM 
and minihalo signals actually cancel each other out if the 
IGM is cold. This is obvious for the bottom curve, but 
it is also true for the others: if minihalos did not exist, 
the rms fluctuations would be larger in both cases. 

6. MINIHALOS AND THE REIONIZATION HISTORY 

Thus, minihalos only dominate the observed fluctua- 
tions if (1) the IGM remains cold and (2) the Lya back- 
ground remains small. How likely are these conditions 
to be fulfilled in realistic structure formation models? In 
this sec tion, we will use the simple models of Furlanetto 
(2006y) to consider minihalos in a global context; we re- 
fer the reader there for a detailed discussion of the vari- 




FlG. 3. — Expected rms brightness temperature fluctuations at 
z = 10 with Tk = T a< j. The solid curve shows the minihalo signal. 
The dashed curves show the net 21 cm fluctuations with J a = 
2 X (K)- 12 , 10 -11 , 10~ 10 ) cm" 2 s- 1 Hz" 1 sr" 1 , from bottom to 
top. 



ous prescriptions. In brief, this model assumes that stars 
dominate the radiation background and sets the total 
star formation rate via the rate at which gas collapses 
onto galaxies, df co \\/dt, where / co ii is the collapse frac- 
tion (or the mass in halos with virial temperatures T vlr > 
10 4 K). The ionizing efficiency is C = ^HeA/escMon, 
where /* is the star formation efficiency, / csc is the es- 
cape fraction of ionizing photons, Ai on is the number of 
ionizing photons produced per baryon incorporated into 
stars, and An e is a normalization constant accounting for 
helium in the IGM. We mod el recombinations following 
IMiralda-Escude et alJ (|2000jh 

The rate of X-ray heating is also assumed to be propor- 
tional to df co \\l dt, calibrated to the local relation between 
star formation rate an d X-ray luminosity, measured be- 
tween 0.2 and 10 key (iGrimm et al~ll2003t iRanalli et all 
I200& IGiifanovetai] I2004J) . Welenore shock heating. 
Finally, we also assume that the soft-UV emissivity (re- 
sponsible for the Wouthuysen-Field effect) is propor- 
tional to df co ii/dt and use the (low-metallicity) Pop- 
ulation II and (high-m ass) Population I I I (henc eforth 
PopII/III)spectral fits of lBarkana & Loebl l|2005bj) to es- 
timate the efficiency with which these photons are pro- 
duced. Note that we properly incorporate h igher L yman- 
series transitions (Hirata 2005; Pritchard & Furlanetto 
2005). 

Although there are obviously a number of free param- 
eters in this model, the general features are easy to un- 
derstand. The 21 cm background contains three major 
transitions: the points when Tk first exceeds T 7 , when 
Lya coupling becomes efficient, and when reionization 
occurs. Because it affects all the backgrounds equally, /+ 
simply shifts everything forward and backward in time. 
X-ray heating only affects the kinetic temperature, while 
/csc only affects the mean ionized fraction x%. The stellar 
initial mass function - especially the choice between Pop 
II and very massive Pop III stellar populations - affects 
the ratio of Lya photons to ionizing photons and moves 
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Fig. 4. — Some characteristics of our example Population II 
history, (a) : The solid line shows Xi , the dotted line shows the col- 
lapsed fraction in star- forming halos / co n, the dashed lines show 
the fraction of gas in minihalos / m h , and the dot-dashed line shows 
bfmh- (b)'. Temperature histories. The thick (or thin) curves in- 
clude (or ignore) X-ray heating. 



Fig. 5. — Minihalo signal in a cosmological context, for a model 
with Pop II stars and without X-ray heating, (a): Fluctuation 
amplitude at k = 0.1 Mpc - 1 . Note that we assume a uniform Xi 
here, (b): Mean brightness temperature relative to the CMB. 



the onset of Lya coupling relative to reionization. Be- 
cause they are so much hotter, Pop III stars push this 
transition later in time, compressing the features in the 
21 cm background. 

We will use a simple Pop II model. We take /* = 
0.1, / osc = 0.1, and N ion = 4000. The solid curve in 
Figure^Ja shows the resulting ionization history. We have 
deliberately chosen /* so that reionization occurs at low 
redshift, thus maximizing the minihalo signal. Note that 
we have ignored feedback processes during reionization; 
they will only push the minihalo era further back in time 
by prolonging reionization. 

Figure 0J> shows the thermal history. The thin solid 
curve in Figure 0J> shows Tk if we ignore X-ray heat- 
ing; the thick solid curve includes it. Note that, even 
with our relatively modest injection rate, X-ray heating 
becomes significant at z ~ 17 and rapidly increases the 
IGM temperature after that point. The dashed curves 
show Ts- Interestingly, Ts — » Tk at z ~ 20: Lya cou- 
pling becomes efficient long before reionization. This is 
because only a small fraction of ionizing photons escape 
their hosts and because Population II stars have a higher 
specific emissivity in the Lyman-line region than in the 
Lyman continuum. 

Figure |ljj also shows the fraction of gas / m h in mini- 
halos. Without heating, / m h increases from ~ 10~ 3 at 
z = 20 to ~ 0.2 at z ~ 6. But with X-ray heating in- 
cluded, it actually peaks at z ~ 15 before plateauing 
somewhat below ~ 0.01. This is, of course, because the 
entropy floor prevents gas from accreting onto the dark 
matter halos. For reference, the dotted curve shows / co n 
(note that we neglect the entropy floor for this quantity) . 

FigurcEJshows the resulting 21 cm histories if we ignore 
X-ray heating. We emphasize that this is not realistic 
and is only meant to build intuition. Panel (a) shows the 



rms temperature fluctuation at k = 0.1 Mpc - , _ while 
panel (b) shows the mean brightness temperature ST),. In 
each case, the dotted, dashed, and solid curves show the 
IGM, minihalo, and total signals. Because the IGM re- 
mains cold but Lya coupling becomes significant at high 
redshifts, we see a net absorption signal that strength- 
ens toward lower redshifts until reionization. Minihalos 
on their own contribute < 0.3 mK in emission; however, 
by removing_gas from the strongly absorbing phase, they 
also reduce 5Tb by a factor / m h- This becomes significant 
at z < 15. 

Because they are highly clustered, minihalos reduce 
the fluctuation amplitude by a larger factor (6/ m h, which 
can be up to ~ 30%; see the dot-dashed curve in Fig. ^ 
which shows this quantity without X-ray heating). Thus, 
in principle, minihalos can be detected indirectly if both 
5Tb and the large-scale rms fluctuation amplitude can 
be measured. If only one of the two measurements is 
available, they are completely degenerate with a differ- 
ent ionized fraction, IGM temperature, or cosmology, 
because the amplitude of each component is quite un- 
certain. These separate measurements, while extremely 
difficult, may in principle be feasible. Radio interferom- 
eters are of course insensitive to the temperature zero 
point, bu t single-dish measu rements may be possible 
(see, e.g., Sha ver et alJll999|) . Also, by careful exami- 
nation of one-point statistics, it ma y be possible to dis - 
cern 5Tb to relatively high precision (Han sen et all2 006) . 
The expected large-scale density fluctuations are easily 
computed from linear theory, provided that one can iso- 
late these fluctuations from any other sources (such as 
ionized bubbles). Alternatively, the fluctuation ampli- 
tude can be constrained via redshift-space distortions 
(Bark ana fc Loeal2005a|) . although in practice that is 

7 This scale is arbitrary and is chosen only because it sits in the 
most easily observed window. The ratio between the curves, which 
is what is relevant to our argument, is independent of scale so long 
as the one-halo term can be ignored. 
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Fig. 6. — Same as Fig.0 except including X-ray heating. The 
dashed curve shows the signal from minihalos with maximal feed- 
back, while the dot-dashed curve shows the signal neglecting their 
suppression. The two solid curves show the net signals in these two 
cases. Note that the IGM and total curves overlap nearly every- 
where, showing explicitly that minihalos are nearly impossible to 
isolate without resolving them. The discontinuity in A corresponds 
to the crossover from absorption to emission. 



rather difficult (McQuinn et al.H2005|) . 

Of course, Figure is far from realistic because we 
must account for X-ray heating. Figure [fj] shows the re- 
sults in such a model. Because the effects of feedback 
remain controversial, we consider two different prescrip- 
tions for it: one in which feedback maximally suppresses 
the minihalo abundance (see |21 dashed lines) and one in 
which it has no effect on them (dot-dashed lines). The 
overall histories are rather different from Fig. with an 
early absorption phase followed by emission once X-ray 
heating kicks in; at the crossover point, both the mean 
signal and the fluctuations nearly vanish in our model 
(note that this is only because we neglect spin tempera- 
ture fluctuations). But most importantly, minihalos are 
now completely swamped by the IGM signal. Lya cou- 
pling turns on at z ~ 17, when / m h ~ 10~ 2 . Beyond 
this point, minihalos become degenerate with the IGM 
fluctuations. In principle, combined measurements of A 
and 5Tb could work when the IGM is in absorption (as 
in Fig. [5J, but by z ~ 15 the entire IGM is hot and 
emits relative to the CMB. Past this point, there is no 
way to separate minihalos from the IGM without resolv- 
ing them because both appear in emission. Interestingly, 
the treatment of the entropy floor makes essentially no 
difference to the results, because feedback kicks in only 
when the IGM is already strongly coupled (and usually 
saturated in emission). 

Thus minihalos arc unlikely to have any observable ef- 
fect on the 21 cm signal (except on small spatial scales, 
obviously). Of course, one can quibble with the particu- 
lar parameters of our model. But we have endeavored to 
choose values that maximize the importance of miniha- 
los, by for example delaying reionization to z ~ 7 and ig- 
noring other sources of fluctuations in ST},. Only by mov- 
ing Lya coupling closer to reionization would matters be 



improved. Perhaps the best possibility is if very massive 
Population III stars dominate the radiation background 
throughout reionization. These have much harder spec- 
tra than normal st ars and so induce later Lya coupling 
(Furlanctto 2006b). But we find that, even in this case, 
minihalos make only a small difference to the signal, and 
then only at z > 12. Moreover, it seems implausible 
that star formation in pristine gas can dominate the ion- 
izing photon budg et throughout reionization (see, e.g., 
iFurlanetto fc Loebll2005]) . 

7. DISCUSSION 

We have presented a new method for calculating the 
21 cm emission fluctuations from minihalos (via the halo 
model). This made it particularly easy to estimate the 
effects of IGM heating (or an 'entropy floor;' OH03) and 
allowed us to compare the minihalo and IGM fluctua- 
tion amplitudes in a self-consistent way. Previous cal- 
culations that only included gravitational physics (and 
hence implicitly assumed both radiative heating and Lya 
pumping of the IGM to be negligible) found that the 
21 cm sig nal from minihalos dominates over that of 
the IGM l|Shapiro et al.ll2005|) . In contrast, we found 
that the minihalo signal can only be unambiguously dis- 
tinguished from the IGM in exceptional circumstances. 
First, the Jeans mass must remain sufficiently low that 
small mass minihalos are able to form; once Tk > 200 K, 
the IGM signal becomes comparable to that of miniha- 
los at z = 10. The required temperature is even smaller 
at higher redshifts. Second, Wouthuysen-Field coupling 
must remain negligible. Once the Lya background be- 
comes significant, the diffuse IGM inevitably shines in ei- 
ther absorption or emission; because it contains so much 
more mass than the minihalo phase, it dominates the 
total signal (|Oh fc Macldl200.'^ . 

More importantly, even if the IGM does remain cold to 
relatively low redshifts z < 10 (which is only marginally 
consistent with the WMAP measurement, because heat- 
ing almost certainly precedes reionization), it is almost 
impossible to distinguish the minihalo contribution from 
that of an IGM with a slightly different reionization his- 
tory. Only if minihalos can be resolved (requiring arcsec- 
ond resolution) will they be robustly separable. 8 Other- 
wise, the best way to infer the_minihalo contribution is 
to measure both the absolute STb(z) and the large-scale 
fluctuation amplitude during a phase in which Tk < T y 
(see Fig. |5J. Minihalos affect the mean signal only 
through their mass fraction, but their fluctuations are 
amplified by their bias. This bias cannot be mimicked 
by a different ionization or thermal history and therefore 
presents a unique signature. However, the cold phase al- 
most certainly ends at z ~ 15 when / m h is still extremely 
small, so taking advantage of this separation my well be 
impossible in practice. Once T$ > T 1 , both phases ac- 
quire the same emission characteristics and become indis- 
tinguishable. In most histories, this later phase is by far 
the longest (and of course occurs preferentially at lower 
redshifts, when minihalos are more common). 

Even if minihalos can be detected in this way (or if the 
diffuse IGM somehow maintains T$ = T 7 for an extended 

8 Nonlinear bias could ease the resolution requirem ents some- 
what , if that can be modeled sufficiently precisely Hliev et all 

2003). 



period of time) , we have shown that minihalo s cannot be 
used to co nstrain cosmological parameters l)Iliev et al.l 
120021 I2003J) because of the inevitable astrophysical un- 
certainties in their emission properties. By far the most 
important is the uncertain thermal history (and hence 
Jeans mass), which strongly affects / m h through the en- 
tropy floor. Instead, in the unlikely event that they are 
seen, minihalos will constrain the astrophysical proper- 
ties of the first luminous sources. 

It is still possible to detect minihalos in other ways, 
and we certainly do not wish to downp l ay their sig- 
nifica nce for reionization ijlliev et alJl2005t iCiardi et al.l 
2005). The most promising possibility is the "21 cm for- 
est," in which minihalos appear as weak 21 cm absorp- 
tion features towa r d a b right backgrou nd radio source 
(ICarilli et alJ t2~002: Fu rlanetto k. Loebll2o5l IFurlanettol 
2006a). Because it is relatively easy to achieve high 
spectral resolution, individual minihalos can be observed 
without any degeneracies. One particularly interesting 



application is to search for neutral regions (identified 
through their 21 cm emission, for example) that lack 
minihalo absorption lines. These would presumably be 
hot regions of the IGM and may be relic H II regions 
(OH03). The main problem with this technique is identi- 
fying sufficiently bright background sources; models pre- 
dict that they s hould exist to z > 10 but will become in - 
creasingly rare (jCarilli et alJl2002HHaiman et alJ l20041. 
Another possibility is to search for non-gaussianitics in- 
duced by the nonlinear gravitational clustering of mini- 
halos, although this is likely to be difficult. 
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